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Abstract 

A quantitative application to real supercooled liquids of the mean-field sce- 
nario for the glass transition {Tg) is proposed. This scenario, based on an anal- 
ogy with spin-glass models, suggests a unified picture of the mode-coupling 
dynamical singularity (Tc) and of the entropy crisis at the Kauzmann tem- 
perature (Tx), with Tc > Tg > Tk- Fitting a simple set of mode-coupling 
equations to experimental light-scattering spectra of two fragile liquids and 
deriving the equivalent spin-glass model, we can estimate not only Tc, but 
also the static transition temperature Tg corresponding supposedly to Tk- 
For the models and systems considered here, Tg is always found above T^, in 
the fluid phase. A comparison with recent theoretical calculations shows that 
this overestimation of the ability of a liquid to form a glass seems to be a 
generic feature of the mean-field approach. 

64.70.Pf, 61.20.Lc, 75.10.Nr 



Typeset using REVT[TiX 



*E-mail: krako@lcp.u-psud.fr 
^E-mail: chalba@lcp.u-psud.fr 



1 



Despite considerable experimental and theoretical efforts, understanding the dynamics 
of supercooled liquids and the related phenomenon of glass transition remains a challenging 
problem of classical statistical mechanics 

Recently, Kirkpatrick et al. [||,^ conjectured that generalized spin-glass models with 
discontinuous one-step replica symmetry breaking (IRSB) transitions (like the spin-glass 
with p-spin interactions (p > 2) or the g-state Potts glass (g > 4)) could be relevant models 
for the description of the structural glass transition. Essentially two points substantiate this 
analogy. Firstly, when considered in the mean-field limit where the interactions between 
spins have infinite range, the generalized spin-glasses display a (static) IRSB transition 
to a spin-glass phase at a temperature that is accompanied by the vanishing of the 
configurational entropy similar to the entropy crisis at Tk hypothesized by Kauzmann for 
glassforming liquids Secondly, in the same mean-field limit the study of the Langevin 
dynamics of such models shows that a dynamical transition takes place at a temperature 
greater than T,. Above T^, the time evolution of the spin correlation function is given by a 
non-linear equation of the same type as those occurring in the ideal mode-coupling theory 
(MCT) for the glass transition of simple liquids, where they describe the time evolution of 
density fluctuations 0, and Td coincides with what is called the critical temperature Tc in 
the context of MCT. Below T^, the systems display a nontrivial free-energy landscape that 
in flnite-range models could lead to slow activated dynamics, as described for instance by 
the Adam-Gibbs theory |0. 

Mean-field models and their extensions to finite-range systems could thus provide a 
consistent framework for the study of the liquid-glass transition, a framework in which 
the mode-coupling approach finds a natural place and that catches important aspects of 
the glass phenomenology. This approach has recently been very fruitful: first-principle 
studies based on the rephca method have been proposed for simple liquid models and 
computer simulation studies of the out-of-equilibrium dynamics of simple liquids |10| show 
aging behaviors qualitatively similar to the one displayed by generalized spin glasses | |11| |. 

A question remains, however, elusive. Indeed, it is known that the dynamical transition 
predicted by the ideal MCT is not experimentally observed: it is 'avoided' because of ergod- 
icity restoring processes not accounted by the theory. In general, it is thought to be replaced 
by a smooth crossover regime where the dynamics changes qualitatively 0. Because of the 
similarities between the ideal MCT and the dynamical aspects of the mean-field theory, 
the same breakdown is expected in the latter when more realistic systems with finite-range 
interactions are considered. Its infiuence on the other aspects of the mean-field picture de- 
scribed above is unknown, but the static transition is usually assumed to survive and the 
dynamical freezing is expected to occur only at the static transition [^. This issue clearly 
deserves further investigation. 

In this Letter, as a step in this direction, we propose a quantitative study of the above 
mean-field scenario for real supercooled liquids. By quantitative, we mean that we will 
extract from experimental data values for the two characteristic temperatures introduced 
within the mean-field approach and discuss these values in relation with known proper- 
ties of the investigated systems. At present, no definite method exists for such a work. 
We propose here to start from a phenomenological mode-coupling approach, the so-called 
schematic approach to experimental data |12|, which is found to be particularly well suited 
for investigating the avoidance of the dynamical transition. Then, taking advantage of the 
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coincidence between the ideal mode-coupling equations and the dynamics of some mean-field 
generalized spin-glass models, the schematic calculation is rephrased in terms of an effec- 
tive spin-glass hamiltonian, whose study allows us to determine, for the glassforming liquids 
under investigation, the location of the two transitions predicted by the mean-field theory. 
Schematic models are simple sets of mode-coupling equations , which have proven to 



be useful in testing the MCT on realistic systems. Indeed, these models can be included 
within a fitting procedure of experimental data over a wide time or frequency range, and 
they allow the calculation of effective mode-coupling parameters (Vertices') describing the 
dynamical evolution of a liquid with varying external conditions ||T^,n-TB|. A major interest 



of these models is that they catch the universal features of the mode-coupling equations, in 
particular the asymptotic scaling results valid near the dynamical transition, but are not 
reduced to them. They make possible to overcome the difficulties and uncertainties arising 
from the need for corrections to the asymptotic critical predictions of the theory and to 
take into account the a relaxation as well as the regime of the high-frequency microscopic 
excitations. The price to pay is that these models are somewhat ad hoc and might possibly 
display non-generic features. 

We make use in this Letter of the results obtained from a previous study [1^ of the 



depolarized light scattering spectra of two so-called 'fragile' glassformers, CKN and salol 
[171 . Only the facts relevant to our present calculation are reviewed here and the reader is 



referred to the corresponding paper for details. The basis of our study is the well studied 
Fi2 model p, defined by the following mode-coupling equation for a correlator (po(t): 

Mt) + ^oMt) + ^iMt) + ^ifm^it - r) 0o(r) dr = 0, 

Jo 

with mo(t) = fi0o(^) + V2<t'1{t), Vi,V2 > 0. To fit the experimental light-scattering spectra 
over the full available frequency range, in particular to reproduce the shape of the peak in 
the THz domain and the extra-intensity of the a-peak, we must add a second correlator (pi, 
whose time evolution is given by a similar equation with mi {t) simply chosen proportional to 
mo{t). The first correlator (f)o{t) accounts in an effective way for the slow modes responsible 
for the slowing down of the relaxation, possibly leading to a dynamical transition. The 
second correlator (pi {t) describes the contribution of additional degrees of freedom that are 
important in the light-scattering processes, but whose slow dynamics is dominated by that 
of (f>o(t). The critical slowing down is thus totally driven by the time evolution of (po, 
whose variation with temperature is entirely encoded in the T-dependence of the 'vertices' 
Vi and V2 (the other parameters are taken T-independent). Because (pi plays no role in 
the parametrization of the dynamical evolution, it can be discarded from the building of 
the effective generalized spin-glass model. The latter is obtained remarking that the long- 
time dynamics described by the F12 equation is identical to that of a mean-field spherical 
spin-glass with spins interacting via both 2-body and 3-body random terms, namely 

l<ii<i2<N l<ii<i2<i3<N 

where the Sj's are spherical spin variables and the couplings Jj^jj and Jiiiaig are independent 
gaussian variables with zero means and variances (Jjiia)^ = Vi/N and (^111213)^ = 2v2/N^ 
||18|| . A remarkable property of these systems is that the correlation function of the Hamil- 
tonian, 
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where gs,s' = s ■ s'/N = ^is'^) /N is the overlap between the spin configurations s and 
s', determines completely both their statics and dynamics |]T9|. With parameters Vi and V2 
extracted from a fit of the dynamics, we are thus able to construct an effective generalized 
spin-glass model allowing us to investigate the mean-field scenario for the glassforming salol 
and CKN. 

The phase diagram of the generalized spin-glass model described above is easily computed 
and is plotted on figure |l[ The dynamical-transition line separating the ergodic domain 
in the vicinity of the origin from the non-ergodic one in the vicinity of infinity has two 
branches: {vi = 1, f2 < 1} corresponding to a continuous or type A transition and {vi = 
2 (^2)^''^ — V2, I < V2 < 4} corresponding to a discontinuous or type B one, respectively. 
At the static level, by using the replica trick to average over the quenched disorder [^], the 
model is exactly solved with a IRSB ansatz, which leads to the following free energy density 

Pfiq, x) = - (1 - x)q') - |(1 - (1 - x)q') - ^ ln(l - (1 - x)q) - ^ ln(l - q), 

where q, the mutual overlap between replicas lying in the same cluster, and x, the cluster 
size, are variational parameters (0 < g, x < 1), with respect to which the free energy has 
to be maximized. Here again, the transition line between the replica-symmetric phase at 
small couplings and the IRSB phase at large couplings has two branches: a continuous 
IRSB transition line coincides with the continuous dynamical transition line, whereas a 
discontinuous IRSB transition line is found beyond the dynamical discontinuous transition 
line. It corresponds to the appearance of a non-zero solution for q with x = 1 and its 
equation as a parametric function of q is given by 

2g2-3g-3(l-g) ln(l - g) - + 2(1 - g) ln(l - g) 

^1 = 2 2n A ' ^2 = 3 — r . 

g2(l-g) g^l-?) 

On figure are also reported the effective vertices obtained from the experimental data 
for the two supercooled liquids. They display two different regimes with varying temper- 
ature. At higher temperatures, they vary linearly and show the apparent evolution of the 
liquids toward the dynamical transition expected from MOT. But, above the corresponding 
transition temperature, the vertices behavior changes: the transition is not reached, and the 
vertices follow the dynamical-transition line without crossing it. The first regime can be 
associated with the domain of validity of the ideal MCT, whereas the second one indicates 
the failure of the theory because of the putative onset of activated processes. Indeed, the 
MCT states that the dynamics is governed by vertices which are purely static quantities and 
thus change smoothly with external parameters. Even if in the case of schematic models the 
connection between the effective vertices and static quantities is somewhat obscured, it is 
reasonable to expect smooth variations of the fitted parameters with temperature. Accord- 
ingly, we concentrate in the following on the high temperature regime and its extrapolation 
to lower temperatures. Note that one can avoid the need for extrapolation using a low 
frequency cut-off for the lowest temperatures . 



With this method, we find the dynamical temperatures = = 257 ± 5 K for salol 
and T(i = Tc = 388 ± 5 K for CKN, in good agreement with previously determined values of 
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the mode- coupling transition temperature. By construction, the mode-couphng transitions 
of the hquids and the dynamical transitions of the effective disordered systems are identical. 

We now turn to the static calculation. We are interested in the discontinuous IRSB 
transition temperature, since it is thought to describe the entropy crisis associated with the 
resolution of the Kauzmann paradox. Indeed, as first pointed by Kauzmann, because the 
heat capacity of a supercooled liquid is substantially greater than that of the underlying 
crystalline solid, reasonable extrapolations of the liquid entropy below the glass transition 
seem to cross the entropy of the crystal at a non-zero temperature T^- To solve this paradox, 
following Gibbs and Di Marzio [^, it is sometimes postulated that, below Tg, a second-order 
transition to an ideal glassy state should exist, at which the configurational contribution to 
the entropy of the liquid vanishes. This mechanism is precisely at work in the mean-field 
models at a discontinuous IRSB transition. Between the dynamical and static transitions, 
one finds that the Gibbs measure is dominated by a number of states exponentially large 
in A^, leading to a finite configurational entropy density (defined as the logarithm of that 
number of states divided by A^). At the static transition, this configurational entropy density 
vanishes and stays zero in the low-temperature phase. 

From the above analysis, we find Ts = 242 ± 5 K for salol and Tg = 376 ± 5 K for 
CKN. These values, which are very close to T^, have to be compared with the experimental 
calorimetric glass transition temperatures [Tg = 220 K for salol, Tg = 333 K for CKN) and 
the empirically determined Kauzmann temperatures {T^ = 175 K for salol [^, no value 
for CKN because the crystalline phase is unstable): Ts is found in both cases above Tg, i.e. 
still in the liquid phase ! All these characteristic temperatures are plotted in figure ^ in the 
case of salol: the total configurational entropy decreases by only 10% between T^ and T,. 

We have investigated the effect of minor modifications of the schematic calculation on 
our results, namely changing the expression of the calculated susceptibility as a functional of 
(po and 01 (each can enter linearly or quadratically in the susceptibility expression) and/or 
changing the second memory function to mi{t) = r (f)o{t) 0i(t). The resulting fits are all of 
equally good quality, and the corresponding vertex trajectories, although slightly different, 
agree with the previous values for both the two transition temperatures and the vertices 
at the dynamical transition. As these changes merely affect (pi, this consistency validates 
our assumption of neglecting it for building the effective model. For more drastic changes 
to the model, in which the F12 equation is replaced by another one (for instance, we tried 
the Fi3 and F29 models where mo(t) = fi0o(^) + "^sipoit) and mo(t) = V2(pl(t) + vgcplit), 
respectively), we have not been able to fit satisfactorily the experimental susceptibilities. 
Indeed, the fitted curves failed to reproduce the location of the susceptibility minimum, its 
shape, or else the position of the a-peak. As these features are crucial for the characterization 
of the dynamics within the MCT framework, the results, in qualitative agreement with those 
obtained with the F12 model, do not appear reliable enough for quantitative purpose. The 
origin of this failure is in general unclear, but in some cases can be related to non-generic 
features of a given schematic model (for instance, the F13 model displays an A3 singularity 
very close to the calculated vertex trajectories). This limitation seems to be severe for our 
static calculation, as it is known from the study of mean-field models like spin-glasses with 
p-spin {p > 2) interactions |^ or the Potts glass ||2^, that the larger the asymptotic value 
of the correlation function at the dynamical transition, the larger the ratio T^/Tg, and the 
F12 model only allows for small values of the former. But we stress that, in our model, the 
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temperature enters in a different way, only through the effective vertices whose dependence 
comes out directly from the fitting procedure to dynamic light-scattering susceptibilities. 
There is, thus, no built-in closeness of and Tg in our work. A related remark is the 
independence of our results with respect to changes in the way (f>o{t) enters the calculated 
susceptibility, thereby showing that its infinite time limit is not a sensitive parameter in the 
study. 

Within our phenomenological implementation of the mean-field approach, we find thus 
that this theory seems to overestimate notably the tendency of a supercooled liquid to form 
a glass. This overestimation shows up at two levels (T^ > Tg on one hand, Tg > Tg on the 
other) with clearly different implications. As stated in introduction, this result is expected 
from the dynamical side of the theory, because of its closeness with the ideal mode-coupling 
theory, whose inadequacy at low temperature is well known, and of the need to take into 
account corrections to mean-field in finite dimension 0. From this point of view, this 
conclusion is not new. What is more unexpected is that the obtained static temperature 
values are clearly located in the fluid domain and do not seem to be associated to any change 
of behavior of the studied systems. These results are thus inconsistent with the arguments 
prescribing that, going beyond mean- field, the real dynamical transition should occur at Tg 
for finite dimensional systems [Q. 

Whether the found overestimation has to be ascribed to an inadequacy of our simple 
phenomenological approach or more fundamentally to the theory itself can not be answered 
here a priori. We can nevertheless remark that theoretical studies of simple liquid models 
based on the replica method tend to support our observations and conclusions. Indeed, 
one finds in general that the location of the static transition agrees well with the glass 
transition found by computer simulation, i.e. obtained with large quenching rates and on 
short observation times. This agreement implies thus an overestimation of the ability of 
the liquids to freeze into a glass on a macroscopic time scale, as found in our calculation. 
Concerning the closeness of the dynamical and static transitions, our results can only be 
compared consistently with theoretical approaches allowing one to compute the locations of 
both transitions within the same framework. This is the case in the papers of Ref. @] only, 
in which the transitions are found very close to each other, just as we find here. 

To summarize, we have investigated a potential application of the mean-field scenario for 
the hquid-glass transition to real supercooled liquids by constructing an effective generalized 
spin-glass model whose dynamics reproduces the standard mode-coupling equations used in 
the ideal MCT of glassforming liquids. We have been able to determine from experimental 
data for two fragile glass-formers both the dynamical and static transitions, which are found 
to be rather close and both located in the liquid phase. This result, in qualitative agreement 
with recent theoretical studies, shows that the theory apparently overestimates the ability of 
a supercooled liquid to freeze into a glass, at least in its simple implementation considered 
here. Whether this deficiency could be cured by employing a more sophisticated (but yet 
unknown) version of the mean-field approach to real supercooled liquids or would require 
a non mean-field treatment accounting explicitly for activated dynamics remains an open 
question. 

We are grateful to Prof. W. Gotze and Dr. M. Fuchs for providing us the codes for 
some parts of the calculations. We are also indebted to Prof. H.Z. Cummins and Dr. G. Li 
for the use of their experimental data. Dr. G. Tarjus and Dr. J. Kurchan are particularly 
acknowledged for fruitful discussions. 
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FIGURES 
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FIG. 1. Phase diagram of the generalized spin-glass model (continuous line: continuous dynam- 
ical and IRSB transitions, dashed line: discontinuous d3mamical transition, dotted line: discontin- 
uous IRSB transition) and effective vertices obtained from the MOT fit (open symbols: CKN, filled 
symbols: salol). Inset: evolution of the vertices with temperature. The superimposed straight lines 
are the linear fits used for the low-temperature extrapolations. 
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FIG. 2. Characteristic temperatures calculated within the mean-field scenario for salol, Tc = 
and Ts, represented on a configurational entropy plot (normalized to the entropy of melting AS"™)- 
Entropy data taken from Ref. |^4|. Also shown are the calorimetric glass transition temperature 
Tg, the melting temperature r^, and the extrapolated Kauzmann temperature Tk- The predicted 
static transition Tg occurs well above Tg. 
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